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Abstract 

The charged-particle production ratios p/p, K~/K + , 7r~/7r + , (p + p)/(vr + + 7r~), 
(K + + K~)/(ir + + 7T~) and (p + p)/(K + + K~) are measured with the LHCb 
detector using 0.3 nb^ 1 of pp collisions delivered by the LHC at \fs = 0.9 TeV 
and 1.8 nb -1 at ^/s = 7 TeV. The measurements are performed as a function of 
transverse momentum p^ and pseudorapidity r/. The production ratios are compared 
to the predictions of several Monte Carlo generator settings, none of which are able 
to describe adequately all observables. The ratio p/p is also considered as a function 
of rapidity loss, Ay = ybeam - 2/, and is used to constrain models of baryon transport. 
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1 Introduction 



All underlying interactions responsible for pp collisions at the Large Hadron Collider 
(LHC) and the subsequent hadronisation process can be understood within the context 
of quantum chromodynamics (QCD). In the non-perturbative regime, however, precise 
calculations are difficult to perform and so phenomenological models must be employed. 
Event generators based on these models must be optimised, or 'tuned', to reproduce 
experimental observables. The observables exploited for this purpose include event vari- 
ables, such as particle multiplicities, the kinematical distributions of the inclusive particle 
sample in each event, and the corresponding distributions for individual particle species. 
The generators can then be used in simulation studies when analysing data to search for 
physics beyond the Standard Model. 

The relative proportions of each charged quasi-stable hadron, and the ratio of an- 
tiparticles to particles in a given kinematical region, are important inputs for generator 
tuning. Of these observables, the ratio of antiprotons to protons is of particular inter- 
est. Baryon number conservation requires that the disintegration of the beam particles 
that occurs in high-energy inelastic non-diffractive pp collisions must be balanced by the 
creation of protons or other baryons elsewhere in the event. This topic is known as baryon- 
number transport. Several models exist to describe this transport, but it is not clear which 
mechanisms are most important in driving the phenomenon [THT3]. Pomeron exchange 
is expected to play a significant role, but contributions may exist from other sources, 
for example the Odderon, the existence of which has not yet been established [T5HTB] . 
Experimentally, baryon-number transport can be studied by measuring p/p, the ratio of 
the number of produced antiprotons to protons, as a function of suitable kinematical 
variables. 

In this paper results are presented from the LHCb experiment for the following pro- 
duction ratios: p/p, K~/K + , 7r~/n + , {p + p)/{n + + (K + + K~)/(n + + 7r~) and 
(p + p)/(K + + K~). The first three of these observables are termed the same-particle ra- 
tios and the last three the different-particle ratios. Only prompt particles are considered, 
where a prompt particle is defined to be one that originates from the primary interaction, 
either directly, or through the subsequent decay of a resonance. The ratios are measured 
as a function of transverse momentum p^ and pseudorapidity 77 = — ln(tan#/2), where 9 
is the polar angle with respect to the beam axis. 

Measurements have been performed of the p/p ratio in pp collisions both at the 
LHC [16], and at other facilities [TTH22] . Studies have also been made of the produc- 
tion characteristics of pions, kaons and protons at the LHC at -y/i = 0.9 TeV at mid- 
rapidity [23]. The analysis presented in this paper exploits the unique forward coverage 
of the LHCb spectrometer, and the powerful particle separation capabilities of the ring- 
imaging Cherenkov (RICH) system, to yield results for the production ratios in the range 
2.5 < 77 < 4.5 at both = 0.9 TeV and y/s = 7 TeV. LHCb has previously published 
studies of baryon transport and particle ratios with neutral strange hadrons [24], and 
results for strange baryon observables at the LHC are also available in the midrapidity 
region [2SH2SI- New analyses have also been made public since the submission of this 
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paper [27] . 

The paper is organised as follows. Section [2] introduces the LHCb detector and the 
datasets used. Section [3] describes the selection of the analysis sample, while Sect. 0] 
discusses the calibration of the particle identification performance. The analysis procedure 
is explained in Sect. [5j The assignment of the systematic uncertainties is described in 
Sect. [6]and the results are presented and discussed in Sect. UJ before concluding in Sect. |8j 
Full tables of numerical results may be found in Appendix [A] Throughout, unless specified 
otherwise, particle types are referred to by their name (e.g. proton) when both particles 
and antiparticles are being considered together, and by symbol (e.g. p or p) when it is 
necessary to distinguish between the two. 

2 Data samples and the LHCb detector 

The LHCb experiment is a forward spectrometer at the Large Hadron Collider with a 
pseudorapidity acceptance of approximately 2 < -q < 5. The tracking system begins 
with a silicon strip Vertex Locator (VELO). The VELO consists of 23 sequential stations 
of silicon strip detectors which retract from the beam during injection. A large area 
silicon tracker (TT) follows upstream of a dipole magnet, downstream of which there are 
three tracker stations, each built with a mixture of straw tube and silicon strip detectors. 
The dipole field direction is vertical, and charged tracks reconstructed through the full 
spectrometer are deflected by an integrated B field of around 4 Tm. Hadron identification 
is provided by the RICH system, which consists of two detectors, one upstream of the 
magnet and the other downstream, and is designed to provide particle identification over 
a momentum interval of 2-100 GeV/c. Also present, but not exploited in the current 
analysis, are a calorimeter and muon system. A full description of the LHCb detector 
may be found in [28] . 

The data sample under consideration derives from the early period of the 2010 LHC 
run. Inelastic interactions were triggered by requiring at least one track in either the 
VELO or the tracking stations downstream of the magnet. This trigger was more than 
99% efficient for all offline selected events that contain at least two tracks reconstructed 
through the whole system. Collisions were recorded both at y/s = 0.9 TeV and 7 TeV. 
During 0.9 TeV running, where the beams were wider and the internal crossing-angle of 
the beams within LHCb was larger, detector and machine safety considerations required 
that each VELO half was retracted by 10 mm from the nominal closed position. For 
7 TeV operation the VELO was fully closed. 

The analysis exploits a data sample of around 0.3nb _1 recorded at y/s = 0.9 TeV 
and 1.8 nb" 1 at y/s = 7 TeV. In order to minimise potential detector-related systematic 
biases, the direction of the LHCb dipole field was inverted every 1-2 weeks of data taking. 
At 0.9 TeV the data divide approximately equally between the two polarities, while at 
7 TeV around two-thirds were collected in one configuration. The analysis is performed 
separately for each polarity. 

The beams collided with a crossing angle in the horizontal plane which was set to 
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compensate for the field of the LHCb dipole. This angle was 2.1 mrad in magnitude at 
y/s = 0.9 TeV and 270 (xrad at \/s = 7 TeV. Throughout this analysis momenta and any 
derived quantities are computed in the centre-of-mass frame. 

Monte Carlo simulated events are used to calculate efficiencies and estimate system- 
atic uncertainties. A total of around 140 million events are simulated at 0.9 TeV and 130 
million events at 7 TeV. The pp collisions are generated by Pythia6.4 (29] and the param- 
eters tuned as described in Ref. [3D]- The decays of emerging particles are implemented 
with the EvtGen package [21], with final state radiation described by Photos [32] • The 
resulting particles are transported through LHCb by Geant4 [33l|3l], which models hits 
in the sensitive regions of the detector as well as material interactions as described in 
Ref. [35] . The decay of secondary particles produced in these interactions is controlled by 
Geant4. Additional Pythia6.4 samples with different generator tunes were produced 
in order to provide references with which to compare the results. These were Perugia 0, 
which was tuned on experimental results from SPS, LEP and the Tevatron, and Peru- 
gia NOCR, which includes an extreme model of baryon transport |36j . 

3 Selection of the analysis sample 

The measurement is performed using the analysis sample, the selection of which is de- 
scribed here. Understanding of the particle identification (PID) performance provided by 
the RICH sample is obtained from the calibration sample, which is discussed in Sect. HI 

Events are selected which contain at least one reconstructed primary vertex (PV) 
within 20 cm of the nominal interaction point. The primary vertex finding algorithm 
requires at least three reconstructed tracksjl 

Tracks are only considered that have hits both in the VELO detector and in the track- 
ing stations downstream of the magnet, and for which the track fit yields an acceptable 
X 2 per number of degrees of freedom (ndf). In order to suppress background from decays 
of long-lived particles, or particles produced in secondary interactions, an upper bound 
is placed on the goodness of fit when using the track's impact parameter (IP) to test 
the hypothesis that the track is associated with the PV (xip < 49). To reduce system- 
atic uncertainties in the calculation of the ratio observables, a momentum cut is imposed 
of p > 5 GeV/c, as below this value the cross-section for strong interaction with the 
beampipe and detector elements differs significantly between particle and anti-particle for 
kaons and protons. If a pair of tracks, % and j, are found to have very similar momenta 
(\pi — Pj | / 1 Pi + Pj| < 0.001), then one of the two is rejected at random. This requirement 
is imposed to suppress 'clones', which occur when two tracks are reconstructed from the 
hit points left by a single particle, and eliminates 0(1%) of candidates. 

The analysis is performed in bins of p^ and t]. In p^ three separate regions are 
considered: < 0.8 GeV/c, 0.8 < Pt < 1.2 GeV/c andpr > 1-2 GeV/c. In r] half-integer 

2 The PV requirement can be approximated in Monte Carlo simulation by imposing a filter at generator 
level which demands at least three charged particles with lifetime cr > 10 -9 m, momentum p > 0.3 GeV/c 
and polar angle 15 < 6 < 460 mrad. 
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bins are chosen over the intervals 3.0 < 77 < 4.5 for p-? < 0.8 GeV/c, and 2.5 < rj < 4.5 
for higher p^ values. The rj acceptance is not constant with because the limited size 
of the calibration samples does not allow for the PID performance to be determined with 
adequate precision below 77 = 3 in the lowest pr bin. The bin size is large compared to 
the experimental resolution and hence bin-to-bin migration effects are negligible in the 
analysis. 

The RICH is used to select the analysis sample at both energy points from which 
the ratio observables are determined. A pattern recognition and particle identification 
algorithm uses information from the RICH and tracking detectors to construct a negative 
log likelihood for each particle hypothesis (e, /i, it, A or p). This negative log likelihood 
is minimised for the event as a whole. After minimisation, the change in log likelihood 
(DLL) is recorded for each track when the particle type is switched from that of the 
preferred assignment to another hypothesis. Using this information the separation in 
log likelihood DLL(x — y) can be calculated for any two particle hypotheses x and y, 
where a positive value indicates that x is the favoured option. In the analysis, cuts are 
placed on DLL(p — A) versus DLL(p — it) to select protons and on DLL (A' — p) versus 
DLL(A' — it) to select kaons. Pions are selected with a simple cut on DLL(tt — A'). As the 
RICH performance varies with momentum and track density, different cuts are applied 
in each (pt,v) bin. The selection cuts are chosen in order to optimise purity, together 
with the requirement that the identification efficiency be at least 10%. Figure [1] shows the 
background-subtracted two-dimensional distribution of DLL(p — A) and DLLQo — tt) for 
protons, kaons and pions in the calibration sample for one example bin. The approximate 
number of positive and negative tracks selected in each PID category is given in Tables [T] 
and 121 A charge asymmetry can be observed in many bins, most noticeably for the 
protons. 



4 Calibration of particle identification 

The calibration sample consists of the decays^) Ag — > tt + tt~, A — > pir" and <\> — > K + K~ , 
all selected from the 7 TeV data. The signal yields in each category are 4.7 million, 1.4 
million and 5.5 million, respectively. 

The Ag and A (collectively termed V^ ) decays are reconstructed through a selection 
algorithm devoid of RICH PID requirements, identical to that used in Ref . [21] , providing 
samples of pions and protons which are unbiased for PID studies. The purity of the 
samples varies across the pi and 77 bins, but is found always to be in excess of 83% and 
87%, for Ag and A, respectively. Isolating <fi — > A + A~ decays with adequate purity is 
only achievable by exploiting RICH information. A PID requirement of DLL(A' — it) > 15 
is placed on one of the two kaon candidates, chosen at random, so as to leave the other 
candidate unbiased for calibration studies. The purity of this selection ranges from 17% 
to 68%, over the kinematic range. Examples of the invariant mass distributions obtained 
in a typical analysis bin for each of the three calibration modes are shown in Fig. [2j 

3 In this section the inclusion of the charge conjugate decay A — > pir + is implicit. 
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Figure 1: Two-dimensional distribution of the change in log likelihood DLL(p— K) and DLL(p— 
7r) for (a) protons, (b) kaons and (c) pions (here shown for negative tracks and one magnet 
polarity) in the calibration sample with px > 1-2 GeV/c and 3.5 < ij < 4.0. The region 
indicated by the dotted lines in the top right corner of each plot is that which is selected in the 
analysis to isolate the proton sample. The selection of the calibration sample is discussed in 
Sect.H 
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Table 1: Number of particle candidates in the analysis sample at y/s = 0.9 TeV, separated into 
positive and negative charge (Q). 







Vt 


< 0.8 GeV/c 


0.8 < p T < 1.2 GeV/c 


T>T > 
/ 1 — 


1.2 GeV/c 




Q 


V 


K 


71 


P 


K 


71 


p 


K 


71 


2.5 < 7] < 3.0 


+ 


— 






16k 


39k 


270k 


19k 


36k 


130k 












13k 


35k 


270k 


13k 


31k 


120k 


3.0 < rj < 3.5 


+ 


21k 


78k 


1.1M 


30k 


63k 


260k 


34k 


39k 


120k 






17k 


69k 


1.1M 


21k 


55k 


250k 


20k 


31k 


100k 


3.5 < T) < 4.0 


+ 


55k 


120k 


1.9M 


55k 


60k 


240k 


31k 


33k 


97k 






38k 


100k 


1.9M 


33k 


49k 


230k 


14k 


23k 


85k 


4.0 < rj < 4.5 


+ 


26k 


90k 


1.2M 


23k 


30k 


100k 


14k 


Ilk 


39k 






21k 


86k 


1.2M 


Ilk 


22k 


88k 


4.2k 


6.6k 


30k 



Table 2: Number of particle candidates in the analysis sample at y/s = 7.0 TeV, separated into 
positive and negative charge (Q) 







Pt 


< 0.8 GeV/c 


0.8 <p T < 1.2 GeV/c 


Pt 


> 1.2 GeV/c 




Q 


V 


K 


71 


P 


K 


71 


P 


K 


71 


2.5 <rj< 3.0 


+ 








59k 


250k 


2.0M 


140k 


360k 


1.3M 












52k 


240k 


2.0M 


130k 


350k 


1.3M 


3.0 < rj < 3.5 


+ 


76k 


451k 


6.6M 


120k 


460k 


1.9M 


240k 


400k 


1.2M 






67k 


420k 


6.6M 


110k 


440k 


1.9M 


210k 


380k 


1.2M 


3.5 < r) < 4.0 


+ 


230k 


730k 


11M 


280k 


450k 


1.8M 


250k 


350k 


1.0M 






200k 


700k 


11M 


240k 


420k 


1.8M 


200k 


320k 


1.0M 


4.0 < 7] < 4.5 


+ 


140k 


950k 


12M 


140k 


370k 


1.3M 


140k 


170k 


740k 






120k 


900k 


12M 


120k 


330k 


1.2M 


110k 


170k 


650k 
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Figure 2: Invariant mass distributions reconstructed for one magnet polarity from the y/s = 
7TeV data in the analysis bin for which the positive final-state particle has > 1-2 GeV and 
3.5 < 7] < 4.0 for (a) K% — > 7r + 7r , (b) A — > pir and (c) <p — > K + K . The results of unbinned 
maximum likelihood fits to the data are superimposed. 



In order to study the PID performance on the unbiased K tracks associated with 
genuine (f> decays the sPlot [37] technique is employed, using the invariant mass as the 
uncorrelated discriminating variable, to produce distributions of quantities such as the 
RICH DLL(K — 7r). Although the background contamination in the V° selections is small 
in comparison, the same strategy is employed to extract the true DLL distributions from 
all unbiased track samples in each analysis bin. The two V° signal peaks are parameterised 
by a double Gaussian function, while the strongly decaying is described by a Breit- 
Wigner function convoluted with a Gaussian. The background is modelled by a first and 
third order Chebyshev polynomial for the V° and <fi distributions, respectively. 

The resulting distributions cannot be applied directly to the analysis sample for two 
reasons. The first is that the PID performance varies with momentum, and the finite 
size of the (pr> v) bins means that the momentum spectrum within each bin is in general 
different between the calibration and analysis samples. The second is that the PID perfor- 
mance is also dependent on multiplicity, and here significant differences exist between the 
calibration and analysis samples, most noticeably for the 0.9 TeV data. To obtain rates 
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applicable to the 0.9 TeV and 7 TeV analysis samples, it is therefore necessary to reweight 
the calibration tracks such that their distributions in momentum and track multiplicity 
match those of a suitable reference sample. A single reference sample cannot be adopted 
for all particle types, as the unbiased momentum spectrum is in general different particle- 
to-particle. For this reason, the analysis samples are used, but with the final selection 
replaced by looser PID requirements. This modified selection minimises distortions to the 
momentum spectra, while providing sufficient purity for the differences in distributions 
between particle species to be still evident. In each (pp; v) bin the reference and cali- 
bration samples are subdivided into six momentum and four track multiplicity cells, and 
the relative proportion of tracks within each cell is used to calculate a weight. The PID 
performance as determined from the calibration samples after reweighting is then applied 
in the analysis. 
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Figure 3: Monte Carlo PID efficiency study for protons (a), kaons (b) and pions (c). Shown is 
a comparison of measured efficiencies from a Monte Carlo calibration sample, after background 
subtraction and reweighting, with the true values in the Monte Carlo analysis sample. The 
diagonal line on each plot is drawn with unit gradient. 



The reliability of the calibration can be assessed by comparing the results for the mea- 
sured PID efficiencies from a Monte Carlo simulated calibration sample, after background 
subtraction and reweighting, to the true values in the Monte Carlo analysis sample. The 
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results are shown in Fig. |3l where each entry comes from a separate (pp> v) bin. In general 
good agreement is observed over a wide range of working points, with some residual biases 
seen at low p^. These biases can be attributed to minor deficiencies in the reweighting 
procedure, which are expected to be most prevalent in this region. 

5 Analysis procedure 

The number of particles, Nf, selected in each of the three classes % — p, K or n, is related 
to the true number of particles before particle identification, N?, by the relationship 



where the matrix element e^j is the probability of identifying particle type i as type j. 
This expression is valid for the purposes of the measurement since the fraction of other 
particle types, in particular electrons and muons, contaminating the selected sample is 
negligible. As Nf and e^j are known, the expression can be inverted to determine A^ T . 
This is done for each (pT,t]) bin, at each energy point and magnet polarity setting. After 
this step (and including the low px scaling factor correction discussed below) the purities 
of each sample can be calculated. Averaged over the analysis bins the purities at 0.9 TeV 
(7 TeV) are found to be 0.90 (0.84), 0.89 (0.87) and 0.98 (0.97) for the protons, kaons 
and pions, respectively. 

In order to relate A^ T to the number of particles produced in the primary interac- 
tion it is necessary to correct for the effects of non-prompt contamination, geometrical 
acceptance losses and track finding inefficiency The non-prompt correction, according 
to simulation, is typically 1-2%, and is similar for positive and negative particles. The 
most important correction when calculating the particle ratios is that related to the track 
finding inefficiency, as different interaction cross-sections and decays in flight mean that 
this effect does not in general cancel. All correction factors are taken from simulation, and 
are applied bin-by-bin, after which the particle ratios are determined. The corrections 
typically lead to a change of less than a relative 10% on the ratios. 

The analysis procedure is validated on simulated events in which the measured ratios 
are compared with those expected from generator level. A \ 2 is formed over all the r] bins 
at low pt, summed over the different-particle ratios. Good agreement is found for the 
same-particle ratios over all rj and px, and for the different-particle ratios at mid and high 
px- Discrepancies are however observed at low px for the different-particle ratios, which 
are attributed to imperfections in the PID reweighting procedure for this region. The 
X 2 in the low px bin is then minimised by applying charge-independent scaling factors of 
1.33 (1.10) and 0.90 (0.86) for the proton and kaon efficiencies, respectively, at 0.9 TeV 
(7 TeV). An uncertainty of ±0.11 is assigned to the scaling factors, uncorrelated bin-to-bin, 
in order to obtain x 2 /ndf 1 at both energy points. This uncertainty is fully correlated 
between positive and negative tracks. Although no bias is observed at mid and high px, an 




(1) 
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additional relative uncertainty of ±0.03 is assigned to the proton and kaon efficiencies for 
these bins to yield an acceptable scatter (i.e. x 2 /ndf ~ 1)- This uncertainty is also taken 
to be uncorrelated bin-to-bin, but fully correlated between positive and negative tracks. 
The scaling factors and uncertainties from these studies are adopted for the analysis of 
the data. 

6 Systematic uncertainties 

The contribution to the systematic uncertainty of all effects considered is summarised in 
Tables [3] and H] for the same-particle ratios, and in Tables [5] and [6] for the different-particle 
ratios. 

The dominant uncertainty is associated with the understanding of the PID perfor- 
mance. Each element in the identification matrix (Eq. [I]), is smeared by a Gaussian 
of width corresponding to the uncertainty in the identification (or misidentification) effi- 
ciency of that element, and the full set of particle ratios is recalculated. This uncertainty is 
the sum in quadrature of the statistical error from the calibration sample after reweight- 
ing, as discussed in Sect. HI and the additional uncertainty assigned after the analysis 
validation, described in Sect. The procedure is repeated many times and the width of 
the resulting distributions is assigned as the systematic uncertainty. As can be seen in 
Tables [3HH1 there is a large range in the magnitude of this contribution. The uncertainty 
is smallest at high p^ and 77, on account of the distribution of the events in the calibration 
sample. For each observable the largest value is found in the lowest r\ bin at mid-pr- If 
this bin and the lowest 77 bin at low p^ are discounted, the variation in uncertainty of the 
remainder of the acceptance is much smaller, being typically a factor of two or three. 

Knowledge of the interaction cross-sections and the amount of material encountered by 
particles in traversing the spectrometer is necessary to determine the fraction of particles 
that cannot be reconstructed due to having undergone a strong interaction. The interac- 
tion cross-sections as implemented in the LHCb simulation agree with measurements [5%] 
over the momentum range of interest to a precision of around 20% for protons and kaons, 
and 10% for pions. The material description up to and including the tracking detectors is 
correct within a tolerance of 10%. The effect of these uncertainties is propagated through 
in the calculation of the track loss for each particle type from strong interaction effects. 

The detection efficiency of positive and negative tracks need not be identical due to the 
fact that each category is swept by the dipole field, on average, to different regions of the 
spectrometer. Studies using J/ijj — > decays in which one track is selected by muon 

chamber information alone constrain any charge asymmetry in the track reconstruction 
efficiency to be less than 1.0 (0.5)% for the 0.9 (7) TeV data. These values are used to 
assign systematic uncertainties on the particle ratios. The identification efficiencies in the 
RICH system are measured separately for each charge, and so this effect is accounted for 
in the inputs to the analysis. A cross-check that there are no significant reconstruction 
asymmetries left unaccounted for is provided by a comparison of the results obtained 
with the two polarity settings of the dipole magnet. Consistent results are found for all 
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observables. 

A possible source of bias arises from the contribution of 'ghost' tracks; these are tracks 
which have no correspondence with the trajectory of any charged particle in the event, 
but are reconstructed from the incorrect association of hit points in the tracking detectors. 
Systematic uncertainties are therefore assigned in each (p?, v) bin for each category of ratio 
by subtracting the estimated contribution of ghost tracks for each particle assignment, and 
determining the resulting shifts in the calculated ratios. A sample enriched in ghost tracks 
can be obtained by selecting tracks where the number of hits associated with the track in 
the TT detector is significantly less than that expected for a particle with that trajectory. 
Comparison of the fraction of tracks of this nature in data and simulation is used to 
determine the ghost-track rate in data by scaling the known rate in simulation. This 
exercise is performed independently for identified tracks which are above and below the 
Cherenkov threshold in the RICH system. The contamination from ghost tracks is lower 
in the above-threshold category since the presence of photodetector hits is indicative of a 
genuine track. The total ghost-track fraction for pions and kaons is found to be typically 
below 1%, rising to around 2% in certain bins. The ghost-track fraction for protons rises 
to 5% in some bins, on account of the larger fraction of this particle type lying below 
the RICH threshold. The charge asymmetry for this background is found to be small 
and the assigned systematic uncertainty is in general around 0.1%. To provide further 
confirmation that ghost tracks are not a significant source of bias the analysis is repeated 
with different cut values on the track-fit xV n df and stable results are found. 

Clones are suppressed by the requirement that only one track is retained from pairs of 
tracks that have very similar momentum. The analysis is repeated with the requirement 
removed, and negligible changes are seen for all observables. 

Contamination from non-prompt particles induces a small uncertainty in the mea- 
surement, as this source of background is at a low level and cancels to first order in the 
ratios. The error is assigned by repeating the analysis and doubling the assumed charge 
asymmetry of these tracks compared with the value found from the simulation. No sig- 
nificant variations are observed when the analysis is repeated with different cut values on 
the prompt-track selection variable Xip- 

The total systematic uncertainty for each observable is obtained by summing in 
quadrature the individual contributions in each (px, v) bin. In general, the system- 
atic uncertainty is significantly larger than the statistical uncertainty, with the largest 
contribution coming from the knowledge of the PID performance, which is limited by the 
size of the calibration sample. 

7 Results 

The measurements of the same-particle ratios are plotted in Figs. HI and El and those 
of the different-particle ratios in Figs. and [9j The numerical values can be found in 
Appendix |A] Also shown are the predictions of several Pythia6.4 generator settings, or 
'tunes': LHCb MC [30], Perugia and Perugia NOCR [36]. At 0.9 TeV the p/p ratio 



11 



Table 3: Range of systematic uncertainties, in percent, for same-particle ratios at y/s = 0.9 
TeV. 





P/p 


K-/K+ 


71 


7tt+ 


PID 


7.5 - 46.7 


4.9 - 42.4 


0.8 


- 6.0 


Cross-sections 


0.2 - 1.6 


0.1 - 1.5 


<0.1 


- 0.8 


Detector material 


0.1 - 0.8 


0.1 - 0.7 


<0.1 


- 0.8 


Ghosts 


<0.1 - 0.1 


<0.1 - 0.1 


<0.1 


- 0.1 


Tracking asymmetry 


1.0 


1.0 




1.0 


Non-prompt 


<0.1 - 0.2 


<0.1 - 0.1 


<0.1 


- 0.1 


Total 


7.7 - 46.7 


5.0 - 42.4 


1.3 


- 6.0 



Table 4: Range of systematic uncertainties, in percent, for same-particle ratios at i/s = 7 TeV. 





p/p 


K-/K+ 


71 


7tt+ 


PID 


3.4 - 26.4 


2.0 - 15.8 


0.6 


- 2.7 


Cross-sections 


0.3 - 1.8 


0.3 - 0.7 


<0.1 


- 0.2 


Detector material 


0.2 - 0.9 


0.1 - 0.4 


<0.1 


- 0.2 


Ghosts 


<0.1 - 0.4 


<0.1 - 0.1 




<0.1 


Tracking asymmetry 


0.5 


0.5 




0.5 


Non-prompt 


<0.1 - 0.2 


<0.1 - 0.1 


<0.1 


- 0.1 


Total 


3.5 - 26.5 


2.1 - 15.8 


0.8 


- 2.8 



Table 5: Range of systematic uncertainties, in percent, for different-particle ratios at y/s = 0.9 
TeV. 





(p+p)/(TT + + 7T ) 


(K+ + K-)/(tt + + 7T-) 


(p + p)/(K+ + K~) 


PID 


10.2 - 63.7 


8.1 - 46.8 


5.9 - 42.6 


Cross-sections 


0.1 - 1.6 


0.4 - 1.3 


0.2 - 2.4 


Detector material 


<0.1 - 0.8 


0.2 - 0.7 


0.1 - 1.2 


Ghosts 


<0.1 - 0.1 


<0.1 - 0.1 


<0.1 - 0.1 


Tracking asymmetry 


<0.1 


<0.1 


<0.1 


Non-prompt 


<0.1 - 0.2 


0.1 


<0.1 - 0.1 


Total 


10.2 - 63.7 


8.6 - 46.8 


6.0 - 42.6 
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Table 6: Range of systematic uncertainties, in percent, for different-particle ratios at y/s = 7 
TeV. 





[p + P)/ [TT + 7T ) 


[K + K ) / {7T + 7T J 


[p + p) 1 + K ) 


PID 


5.9 - 31.1 


4.6 - 26.6 


3.7 - 16.1 


Cross-sections 


0.3 - 2.2 


1.2 - 1.5 


0.2 - 2.1 


Detector material 


0.2 - 1.1 


0.6 - 0.8 


0.1 - 1.0 


Ghosts 


<0.1 - 0.3 


<0.1 - 0.3 


<0.1 - 0.2 


Tracking asymmetry 


<0.1 


<0.1 


<0.1 


Non-prompt 


<0.1 - 0.3 


0.1 - 0.2 


<0.1 - 0.2 


Total 


6.0 - 31.1 


4.8 - 26.7 


3.7 - 16.2 



falls from around 0.8 at low 77 to around 0.4 in the highest pt and 77 bin. At this energy 
point there is a significant spread between models for the Monte Carlo predictions, with 
the data lying significantly below the LHCb MC and Perugia expectations, but close to 
those of Perugia NOCR. At higher energy the p/p ratio is higher and varies more slowly, 
in good agreement with LHCb MC and Perugia and less so with Perugia NOCR. The 
K~ /K + and 7r~ /tt + ratios also differ from unity, most noticeably at high p^ and high 77. 
This behaviour is in general well modelled by all the generator tunes, which give similar 
predictions for these observables. Small discrepancies are observed at 7 TeV for K~ / K + 
at low pt, and tt~ /tt + at high p^. When comparing the measurements and predictions 
for the different-particle ratios the most striking differences occur for (p + p)/{tt + + 7r~) 
and (K + + K~)/{ti + + 7r~), where there is a tendency for the data to lie significantly 
higher than the Perugia and NOCR expectations. The agreement with the LHCb MC 
for these observables is generally good. 

It is instructive to consider the p/p results as a function of rapidity loss, Ay = ybca,m—y, 
where ?/bcam is the rapidity of the protons in the LHC beam which travels forward in 
the spectrometer (|/beam = 6.87 at 0.9 TeV and 8.92 at 7 TeV). For the same-particle 
ratios it is possible to determine the rapidity value to which the measurement in each 
77 bin corresponds. In each bin the mean and RMS spread of the rapidity of the tracks 
in the analysis sample is determined. Correlations are accounted for, but these are in 
general negligible as the uncertainties are dominated by the PID errors, which for these 
observables are statistical in nature. A small correction is applied to this mean, obtained 
from Monte Carlo, to account for the distortion to the unbiased spectrum that is induced 
by the reconstruction and PID requirements. The values of the mean and RMS spread 
of the rapidities for p/p can be found in Appendix [Al together with those of K~/K + 
and 7T _ /tt + . As no evidence is seen of any p^ dependence in the distribution of the p/p 
results against Ay the measurements in each 77 bin at each energy point are integrated 
over p^, with the uncertainties on the individual values of the ratios used to determine the 
weights of each input entering into the mean. The mean p/p ratios are given as a function 
of Ay in Table [7] and plotted in Fig. [10], with the results from other experiments [T6T - 
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Table 7: Results for p/p ratio integrated over pt in V bins as a function of the rapidity loss Ay. 





r\ range 


Ay 


Ratio 


0.9 TeV 


4.0-4.5 


3.1 ±0.2 


0.48 ±0.03 




3.5-4.0 


3.5 ±0.2 


0.57 ±0.02 




3.0-3.5 


3.9 ±0.2 


0.65 ±0.03 




2.5-3.0 


4.3 ±0.1 


0.81 ±0.09 


7TeV 


4.0-4.5 


5.1 ±0.2 


0.90 ±0.03 




3.5-4.0 


5.5 ±0.2 


0.92 ±0.02 




3.0-3.5 


5.9 ±0.2 


0.91 ±0.02 




2.5-3.0 


6.3 ±0.1 


0.89 ±0.04 



121] superimposed. The LHCb results cover a wider range of Ay than any other single 
experiment and significantly improve the precision of the measurements in the region 
Ay < 6.5. 

Within the Regge model, baryon production at high energy is driven by 
Pomeron exchange and baryon transport by string-junction exchange [9]. Assum- 
ing this picture the Ay dependence of the p/p ratio approximately follows the form 
1/ (1 ± Cexp[(«j — ap)Ay]), where C determines the relative contributions of the two 
mechanisms, and aj (ap) is the intercept of the string junction (Pomeron) Regge tra- 
jectory. Figure [TU] shows the results of fitting this expression to both the LHCb and, in 
order to constrain the high Ay region, the ALICE data. Both C and (aj — ap) are free 
parameters of the fit and are determined to be 22.5 ± 6.0 and —0.98 ± 0.07 respectively 
with a x 2 / n df of 8.7/8. Taking ap = 1.2 [39] suggests a low value of aj, significantly 
below the a,j ~ 0.5 expected if the string-junction intercept is associated with that of 
the standard Reggeon (or meson). The value of aj ~ 0.9 which would be expected if 
the string junction is associated with the Odderon [T5] is excluded using this fit model. 
The same conclusion applies if the LHCb and ALICE p/p ratio values are fitted with an 
alternative parameterisation [TT] C ■ (s[GeV 2 ]Y aj ~ ap ^ 2 • cosh[y(aj — ap)}, which yields 
the results C = 10.2 ± 1.8, (atj - a P ) = -0.86 ± 0.05 with a x 2 /ndf of 10.2/8. 

8 Conclusions 

Measurements have been presented of the charged-particle production ratios p/p, K~ /K + , 
ix~ /tt + , (p + p)/(K + + K~), (K + + K~)/(ir + + 7T") and (p + p)/(ir + + ti~) at both 
= 0.9 TeV and ^/s = 7 TeV. The results at 7 TeV are the first studies of pion, kaon 
and proton production to be performed at this energy. Comparisons have been made 
with several generator tunes (LHCb MC, Perugia and Perugia NOCR). No single tune 
is able to describe well all observables. The most significant discrepancies occur for the 
(p + p)/ (7r + + 7r~) and (K + ± K~)/ [ti + + 7r _ ) ratios, where the measurements are much 
higher than the Perugia and Perugia NOCR predictions, but lie reasonably close to the 
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Figure 4: Results for the p/p ratio at 0.9 TeV (a) and 7 TeV (b). 
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Figure 5: Results for the K~/K+ ratio at 0.9 TeV (a) and 7 TeV (b). 
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Figure 6: Results for the 7r-/vr + ratio at 0.9 TeV (a) and 7 TeV (b). 
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Figure 7: Results for the (p + p)/(vr+ + vr") ratio at 0.9 TeV (a) and 7 TeV (b). 
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Figure 8: Results for the (K+ + K~)/(tt+ + vr~) ratio at 0.9 TeV (a) and 7 TeV (b). 
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Figure 9: Results for the (p + p)/(K+ + K~) ratio at 0.9 TeV (a) and 7 TeV (b). 
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Figure 10: Results for the p/p ratio against the rapidity loss Ay from LHCb. Results from 
other experiments are also shown [16H21] , Superimposed is a fit to the LHCb and ALICE [16| 
measurements that is described in the text. 
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LHCb MC expectation. 

The p/p ratio has been studied as a function of rapidity loss, Ay. The results span the 
Ay interval 3.1 to 6.3, and are more precise than previous measurements in this region. 
Fitting a simple Regge theory inspired model to the LHCb measurements, and those from 
the midrapidity region obtained by ALICE [16], yields a result with a string-junction 
contribution with low intercept value. 

These results, together with those for related observables obtained by LHCb [21], will 
help in understanding the phenomenon of baryon-number transport, and the development 
of hadronisation models to improve the description of Standard Model processes in the 
forward region at the LHC. 
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Appendix 

A Tables of results 

The results for the same-particle ratios, including the rapidity to which the events in 
each pseudorapidity bin correspond, are given in Tables El [9] and [TOj The results for the 
different-particle ratios can be found in Tables [TT] [T2] and [13] 
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Table 8: Results for the p/p ratio with statistical and systematic uncertainties, as a function of pt and r\. Also shown is the mean 
rapidity, y, and RMS spread for the sample in each rj bin. 





p T < 0.8 GeV/c 
y (RMS) Ratio 


0.8 < p T < 1.2 GeV/c 
y (RMS) Ratio 


p T > 1.2 GeV/c 
y (RMS) Ratio 


y/s = 0.9 TeV 


2.58 (0.27) 0.751 ± 0.011 ± 0.163 
2.96 (0.11) 0.729 ± 0.007 ± 0.040 
3.34 (0.24) 0.660 ± 0.009 ± 0.046 


2.42 (0.24) 1.107 ± 0.020 ± 0.349 
2.96 (0.25) 0.684 ± 0.010 ± 0.049 
3.40 (0.22) 0.576 ± 0.007 ± 0.032 
3.87 (0.14) 0.451 ±0.009 ±0.038 


2.63 (0.16) 0.794 ± 0.015 ± 0.089 
3.08 (0.23) 0.614 ± 0.010 ± 0.047 
3.56 (0.24) 0.456 ± 0.009 ± 0.033 
4.02 (0.25) 0.328 ± 0.010 ± 0.049 


2.5 < 7] < 3.0 
3.0 < 7] < 3.5 
3.5 < 7] < 4.0 
4.0 < r] < 4.5 


y/s — 7 TeV 


2.55 (0.27) 0.734 ± 0.011 ± 0.124 
2.96 (0.09) 1.015 ± 0.009 ± 0.037 
3.34 (0.21) 0.957 ± 0.010 ± 0.051 


2.41 (0.25) 1.181 ± 0.020 ± 0.195 
2.98 (0.25) 0.942 ± 0.011 ± 0.036 
3.40 (0.23) 0.916 ± 0.007 ± 0.022 
3.86 (0.19) 0.906 ± 0.010 ± 0.039 


2.63 (0.16) 0.880 ± 0.009 ± 0.039 
3.12 (0.22) 0.905 ± 0.008 ± 0.026 
3.59 (0.24) 0.903 ± 0.008 ± 0.023 
4.06 (0.25) 0.831 ± 0.010 ± 0.050 


2.5 < T) < 3.0 
3.0 < rj < 3.5 
3.5 < r] < 4.0 
4.0 < rj < 4.5 



Table 9: Results for the K /K + ratio with statistical and systematic uncertainties, as a function of and V- Also shown is the 
mean rapidity, y, and RMS spread for the sample in each rj bin. 





p T < 0.8 GeV/c 
?/ (RMS) Ratio 


0.8 < p T < 1.2 GeV/c 
y (RMS) Ratio 


p T > 1.2 GeV/c 
y (RMS) Ratio 


y/s = 0.9 TeV 


2.99 (0.25) 0.834 ± 0.007 ± 0.069 
3.32 (0.25) 1.001 ± 0.007 ± 0.064 
3.67 (0.18) 1.002 ± 0.007 ± 0.093 


2.65 (0.19) 0.870 ± 0.010 ± 0.267 
3.12 (0.21) 0.847 ± 0.009 ± 0.040 
3.62 (0.22) 0.792 ± 0.009 ± 0.028 
4.11 (0.25) 0.680 ± 0.010 ± 0.041 


2.69 (0.14) 0.936 ± 0.013 ± 0.069 
3.18 (0.15) 0.783 ± 0.011 ± 0.037 

3.70 (0.17) 0.723 ± 0.012 ± 0.031 
4.20 (0.21) 0.506 ± 0.014 ± 0.050 


2.5 < r) < 3.0 
3.0 < rj < 3.5 
3.5 < rj < 4.0 
4.0 < 7] < 4.5 


yfs = 7 TeV 


3.02 (0.25) 0.992 ± 0.006 ± 0.063 
3.34 (0.25) 1.062 ± 0.005 ± 0.040 
3.72 (0.22) 1.161 ± 0.005 ± 0.055 


2.65 (0.19) 0.995 ± 0.008 ± 0.101 
3.12 (0.21) 0.966 ± 0.006 ± 0.019 
3.62 (0.21) 0.948 ± 0.006 ± 0.014 
4.11 (0.23) 0.898 ± 0.006 ± 0.025 


2.70 (0.13) 0.991 ± 0.007 ± 0.021 

3.20 (0.14) 0.999 ± 0.006 ± 0.016 
3.70 (0.15) 0.930 ± 0.006 ± 0.017 

4.21 (0.18) 0.958 ± 0.009 ± 0.049 


2.5 < 7] < 3.0 
3.0 < 7} < 3.5 
3.5 < 7] < 4.0 
4.0 < rj < 4.5 



Table 10: Results for the ir /tt + ratio with statistical and systematic uncertainties, as a function of pt and r/. Also shown is the 
mean rapidity, y, and RMS spread for the sample in each rj bin. 





p T < 0.8 GeV/c 
y (RMS) Ratio 


0.8 < p T < 1.2 GeV/c 
y (RMS) Ratio 


p T > 1.2 GeV/c 
y (RMS) Ratio 


y/s = 0.9 TeV 


3.23 (0.09) 0.979 ± 0.005 ± 0.010 
3.71 (0.15) 0.968 ± 0.004 ± 0.011 
4.15 (0.24) 0.929 ± 0.004 ± 0.017 


2.74 (0.07) 0.987 ± 0.010 ± 0.013 
3.23 (0.07) 0.971 ± 0.011 ± 0.010 

3.75 (0.08) 0.951 ± 0.012 ± 0.010 
4.30 (0.10) 0.971 ± 0.016 ± 0.019 


2.75 (0.05) 0.970 ± 0.016 ± 0.014 
3.24 (0.05) 0.926 ± 0.017 ± 0.014 
3.75 (0.05) 0.871 ± 0.019 ± 0.012 
4.30 (0.07) 0.816 ± 0.025 ± 0.029 


2.5 < r) < 3.0 
3.0 < rj < 3.5 
3.5 < rj < 4.0 
4.0 < 7] < 4.5 


yfs = 7 TeV 


3.23 (0.09) 1.011 ± 0.004 ± 0.006 
3.70 (0.14) 1.002 ± 0.003 ± 0.006 
4.14 (0.22) 0.976 ± 0.003 ± 0.006 


2.74 (0.07) 1.002 ± 0.007 ± 0.006 
3.24 (0.07) 0.998 ± 0.007 ± 0.004 
3.74 (0.07) 1.003 ± 0.008 ± 0.004 
4.26 (0.08) 0.998 ± 0.009 ± 0.008 


2.74 (0.04) 1.015 ± 0.010 ± 0.005 
3.24 (0.04) 0.998 ± 0.010 ± 0.004 

3.75 (0.05) 1.000 ± 0.011 ± 0.005 
4.26 (0.05) 0.974 ± 0.012 ± 0.017 


2.5 < 7] < 3.0 
3.0 < 7} < 3.5 
3.5 < 7] < 4.0 
4.0 < rj < 4.5 



Table 11: Results for the {p + p)/{^ + ±7r ) ratio with statistical and systematic uncertainties, 
as a function of px and V- 





p T < 0.8 GeV/c 


0.8 < p T < 1.2 GeV/c 


pr > 1-2 GeV/c 


y/s = 0.9 TeV 


0.086 ±0.001 ±0.021 
0.062 ±0.001 ±0.008 
0.076 ±0.001 ±0.010 


0.328 ±0.007 ±0.104 
0.208 ±0.004 ±0.016 
0.175 ±0.003 ±0.011 
0.233 ±0.006 ±0.022 


0.300 ±0.008 ±0.034 
0.272 ±0.007 ±0.023 
0.252 ±0.007 ±0.020 
0.301 ±0.013 ±0.047 


2.5 < 7] < 3.0 
3.0 < 7] < 3.5 
3.5 < 7] < 4.0 
4.0 < r) < 4.5 


y/s — 7 TeV 


0.085 ±0.001 ±0.017 
0.069 ±0.001 ±0.008 
0.051 ±0.001 ±0.007 


0.235 ±0.004 ±0.039 
0.174 ±0.002 ±0.009 
0.156 ±0.002 ±0.006 
0.184 ±0.003 ±0.010 


0.262 ±0.004 ±0.014 
0.245 ±0.003 ±0.011 
0.242 ±0.003 ±0.010 
0.244 ±0.004 ±0.017 


2.5 < r) < 3.0 
3.0 < r) < 3.5 
3.5 < rj < 4.0 
4.0 < 7] < 4.5 



Table 12: Results for the {K + + K )/(vr + + ir ) ratio with statistical and systematic uncer- 
tainties, as a function of pr and ?]. 





p T < 0.8 GeV/c 


0.8 < p T < 1.2 GeV/c 


p T > 1.2 GeV/c 


^ = 0.9 TeV 


0.180 ±0.002 ±0.026 
0.171 ±0.001 ±0.023 
0.173 ±0.001 ±0.025 


0.184 ±0.003 ±0.056 
0.267 ±0.004 ±0.015 
0.247 ±0.004 ±0.011 
0.268 ±0.006 ±0.018 


0.351 ±0.008 ±0.028 
0.319 ±0.008 ±0.018 
0.314 ±0.009 ±0.017 
0.281 ±0.012 ±0.031 


2.5 < 7] < 3.0 
3.0 < 7] < 3.5 
3.5 < 7] < 4.0 
4.0 < 7] < 4.5 


yfs = 7 TeV 


0.181 ±0.001 ±0.024 
0.173 ±0.001 ±0.021 
0.131 ±0.001 ±0.016 


0.224 ±0.002 ±0.024 
0.263 ±0.003 ±0.010 
0.262 ±0.003 ±0.009 
0.275 ±0.003 ±0.011 


0.371 ±0.004 ±0.014 
0.357 ±0.004 ±0.012 
0.367 ±0.005 ±0.013 
0.328 ±0.005 ±0.020 


2.5 < 77 < 3.0 
3.0 < 7] < 3.5 
3.5 < rj < 4.0 
4.0 < 7] < 4.5 
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Table 13: Results for the (p+p)/(K + + K ) ratio with statistical and systematic uncertainties, 
as a function of p<r and T]. 





p T < 0.8 GeV/c 


0.8 <pr< 1.2 GeV/c 


p T > 1-2 GeV/c 


y/s = 0.9 TeV 


0.481 ±0.008 ±0.139 
0.363 ±0.004 ±0.066 
0.433 ±0.007 ±0.086 


1.831 ±0.039 ±0.822 
0.779 ±0.014 ±0.073 
0.709 ±0.012 ±0.055 
0.865 ±0.021 ±0.097 


0.855 ±0.020 ±0.119 
0.851 ±0.019 ±0.084 
0.799 ±0.021 ±0.076 
1.067 ±0.045 ±0.200 


2.5 < 7] < 3.0 
3.0 < r) < 3.5 
3.5 < r) < 4.0 
4.0 < r) < 4.5 


y/s = 7 TeV 


0.465 ±0.008 ±0.111 
0.398 ±0.004 ±0.067 
0.379 ±0.004 ±0.068 


1.051 ±0.020 ±0.204 
0.660 ±0.009 ±0.039 
0.593 ±0.006 ±0.031 
0.671 ± 0.009 ± 0.046 


0.705 ±0.009 ±0.046 
0.682 ±0.007 ±0.038 
0.659 ±0.007 ±0.037 
0.744 ±0.011 ±0.069 


2.5 < r) < 3.0 
3.0 < 7] < 3.5 
3.5 < rj < 4.0 
4.0 < 7] < 4.5 
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